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Abstract 

More than ten years ago a first step towards quantum error correction (QEC) was implemented 
[Phys. Rev. Lett. 81, 2152 (1998)]. The work showed there was sufficient control in nuclear 
magnetic resonance (NMR) to implement QEC, and demonstrated that the error rate changed from 
e to approximatively ~ e 2 . In the current work we reproduce a similar experiment using control 
techniques that have been since developed, such as GRAPE pulses. We show that the fidelity 
of the QEC gate sequence, and the comparative advantage of QEC are appreciably improved. 
This advantage is maintained despite the errors introduced by the additional operations needed to 
protect the quantum states. 

PACS numbers: 03.67.Lx 
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I. INTRODUCTION 

Quantum computers could solve some problems faster than classical computers [1[. P er ~ 
forming a quantum computation relies on the ability to preserve the coherence of quantum 
states long enough for gates composing the algorithm to be implemented. In practice, the 
quantum coherence is sensitive to the uncontrolled environment and easily damaged by the 
interactions with the environment, a process called decoherence [^J. To protect the fragile 
quantum coherence needed for quantum computation, schemes of quantum error correction 
(QEC) and fault- tolerant quantum computation have been developed Jsj] . 
The 3-bit QEC code was implemented in a liquid-state NMR quantum information processor 
in 1998 as the first experimental demonstration of QEC |4|. More recently, it has been 
implemented in trapped ion and solid-state systems p, @j. Here we report on using the 
GRAPE algorithm |7| to implement a high fidelity version of the 3-bit QEC code for phase 
errors in liquid state NMR. The errors due to natural transversal relaxation are shown to 
be suppressed to a first order. In comparison with the work performed in 1998 the pulse 
sequence fidelity is improved by about 20%, and the reduction of the first order in the decay 
of the remaining polarization after error correction is improved by a factor of ~ 2.3. The 
advantage of the QEC is obtained although the extra operations for protecting the quantum 
states in QEC are subject to errors in implementation. 



II. EXPERIMENTAL PROCEDURE AND RESULTS 

In the current implementation, we use 13 C labelled trichloroethylene (TCE) dissolved in 
d-chloroform as the sample. Data were taken with a Bruker DRX 700 MHz spectrometer. 
The structure of the molecule and the parameters of the spin qubits are shown in Fig. [TJ 
where we denote H as qubit 1, Ci as qubit 2 and C 2 as qubit 3. The Hamiltonian of the 
three-spin system can be written as 
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H = _ 7T J2 ViZi + ^[J 12 ZiZ 2 + J l3 Z,Z 3 + J 23 (X 2 X 3 + Y 2 Y 3 + Z 2 Z 3 )} (1) 
i=i 

where Xi, Yi, Zi denote the Pauli matrices with i indicating the spin location, z/j denotes the 
chemical shift of spin i, and Jij denotes the spin-coupling between spins i and j. The two 
carbon spins are treated in the strongly coupled regime, because the difference in frequencies 
between the two carbons is not large enough for the weak coupling approximation [8|. 
We exploit radio- frequency (r.f.) spin selection techniques to improve the linewidth, and 



9I [lo|. The effect of pulse imperfections due to 



hence the coherence, of the ensemble qubits 
r.f. inhomogeneities is reduced by spatially selecting molecules from a small region in the 
sample through the r.f. power. We choose Ci as the qubit to carry the state for encoding 
and the output state after decoding and error correction. The labelled pseudo-pure states 
0X0 and 0F0, used as the reference states with blank ancilla, are prepared by the circuit in 
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Ref. [10|], where the order is arranged as qubits 1 to 3 and = |0)(0|. The qubit readout is 
performed on Ci, and the signals are normalized with respect to 0X0 or 0Y0, for different 
input states. 

The quantum network used for implementing the QEC code is shown as Fig. [2] (a), where 
Pi n is chosen as X, Y and Z, in separate sequences. We optimize the encoding operation, 
and the decoding operation combined with the error correction as two GRAPE pulses [t| 
with theoretical fidelity > 99.9%. To test the ability of the code to correct for the natural 
dephasing errors due to the transversal relaxation of the spins, the internal spin Hamiltonian 
dTJ is refocused during the time delay implemented between the encoding and decoding 
processes. The refocusing pulse sequence is shown in Fig. [2] (b) where the selective ir pulses 
applied to spin H are hard rectangle pulses with a duration of 20 /is, while the 7r pulses 
applied to Ci or C2 are GRAPE pulses with a duration of 2 ms. Taking into account the 
strong coupling in the Hamiltonian (pQ), we choose the phases of the 7r pulses shown in Fig. 
[2] (b) to obtain a fidelity \Tr{U re f ocus E}\/8 > 99.96%, where U re f ocus denotes the simulated 
unitary implemented by the pulse sequence, and E denotes the identity operation. 
We choose the input states as pi n — X, Y and Z, and measure the polarization that remains 
after error correction in p out . The polarization ratios are denoted as f x , f y and f z . We use 
"entanglement fidelity", represented as 

f = (l + f x + f y + f z )/A (2) 
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to characterize how well the quantum information in p in is preserved 
The experimental results for QEC are shown in Fig. [3] (a). For each delay time, five exper- 
iments are repeated in order to average the random experimental errors in implementation. 
The results of error correction (EC) are represented by •. By averaging the points for each 
delay time, we obtain the averaged entanglement fidelity / shown as x , which can be fitted 
to 0.9828 - 0.0166* - 0.5380* 2 + 0.0014* 3 , with relative fitting error 0.73%, shown as the 
thick dash-dotted curve. 

In order to estimate the performance of the error correction for the encoded states, we 
calculate the entanglement fidelity of decoding (DE) through measuring the remaining po- 
larization before the application of the Toffoli gate, used as the error- correcting step. In this 
case, the decoding operation is implemented by one GRAPE pulse with theoretical fidelity 
> 99.9%. Similar to the measurement for error correction, we also repeat five experiments 
for each delay time. The results are shown as O i n Fig. El (a), and the data points after 
average are marked by +, which can be fitted as 0.9982 - 0.4361* + 0.1679* 2 + 0.2152t 3 , 
with relative fitting error 0.57%, shown as the thick solid curve. Here the ratio of the first 
order decay terms for the two fits is found to be 26.2 ± 0.3. The important reduction of 
the first order decay term indicates the high quality of state stabilization provided by QEC. 
As a comparison, we include the experimental data from Ref. 0], which are marked as 
o and □ in Fig. [3] (b) for the results of QEC and decoding. The data can be fitted as 
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0.7895 - 0.0957* - 0.0828* 2 + 0.0370* 3 and 0.8539 - 1.1021* + 0.8696t 2 + 0.0378t 3 with rela- 
tive fitting errors 0.89% and 0.98%, respectively. The ratio of the first-order decay terms is 
11.5 ±0.2. 

In implementing the QEC code, the operations associated with encoding, decoding and error 
correction are subject to errors, which would lower the ability of the code to protect the 
quantum states. To estimate the effects of the errors, we measure the free evolution decay 
(FED) of pi n under the refocusing sequence shown in Fig. [2] (b). Five experiments are 
repeated for each delay time, and the experimental data for / are shown as A in Fig. [3] (a). 
The averaging points, shown as *, can be fitted as 1.0056 - 0.4164* + 0.3363* 2 - 0.2123t 3 
with relative fitting error 0.45%, shown as the dashed curve. The ratio of the first order 
decay terms in the fits of FED and EC is 25.0 ± 0.3. Through comparing the results of 
QEC and FED, one can find that the errors removed by the QEC code can exceed the errors 
introduced by the extra operations required by the code for delay time > 0.0672 s (~ 6% of 
Ci'sTa). 



III. DISCUSSION 

The pulse durations for encoding, decoding, and the combination of decoding and error 
correction are 8 ms, 8ms, and 13.6 ms, respectively. We exploit the results from simulation 
with ideal pulses to estimate the errors due to the imperfection in pulse implementation. In 



simulation, we choose an uncorrelated error model for T2 errors and ignore T\ errors [12]. 
We represent the measured fidelity as f — Af ideah where f ideal denotes the ideal fidelity 
by simulation and A denotes a factor to estimate the deviation between experiment and 
simulation. One should note that the theoretical entanglement fidelity of DE is the same as 



FED [12j. By fitting the data, we obtain A = 0.983±0.006, 0.998±0.007 and 1.0098±0.0064 
for EC, DE and FED, respectively. The fitting results are shown as the thin dash-dotted, 
solid and dashed curves in Fig. [3] (a). From the simulation results, we estimate the errors 
in implementing the operations associated with the QEC codes are about 1.2% for DE and 
2.7% for EC. 



IV. CONCLUSION 

We optimize the encoding, decoding, and error correction as GRAPE pulses with high the- 
oretical fidelities (> 99.9%). The refocusing sequence is exploited to suspend the evolution 
of the Hamiltonian ([1]) with high fidelity (> 99.96%). The quality of readout signals are 
further improved by r.f. selection. Compared with the experimental results of QEC obtained 
in 1998 [4J], the pulse sequence fidelity is improved by about 20%. By the comparison with 
the free evolution decay, one can benefit from QEC even when errors exist in implementing 
the operations required for QEC. The improvement provided by the error correction is also 
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demonstrated by the reduction of the first order in the decay of the remaining polarization 
after error correction, compared with the decay of the encoded states recovered by decoding 
and free evolution decay of the input states. The "QEC advantage" for the encoded states 
is improved by a factor of ~ 2.3 from the 1998 result. In the current experiment, the second 
order term in the decay after error correction is larger than the previous experiment, be- 
cause of the larger phase errors due to the shorter T 2 time constants [see Fig. Q] (b) , noting 
that T 2 's are 3 s for H, 1.1 s for Ci and 0.6 s for C 2 in the previous experiment 4j. The 
experimental errors arise mainly from the imperfection in implementing the GRAPE pulses. 
Additionally inhomogeneities of magnetic fields and the limitation of Ti also contribute to 
errors. 
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FIG. 1: Parameters of the spin qubits. (a) Chemical shifts shown as the diagonal terms and the 
couplings between spins shown as the non-diagonal terms in Hz. The inset shows the molecule 
structure where the three qubits are H, Ci and C2. (b) The relaxation times Ti's are measured by 
the standard inversion recovery sequence. T2's are measured by the Hahn-echo with one refocusing 
pulse, by ignoring the strong coupling in the Hamiltonian ([T]). 
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FIG. 2: (Color online). Quantum network for quantum error correction (a) where the T2 noise 
is introduced by a variable time delay implemented by the pulse sequence (b) which refocuses 
the evolution of the Hamiltonian (JTJ) to an identity operation with theoretical fidelity higher than 
99.96%. In (b) the unfilled rectangle represents a hard tt pulse with duration of 20 /is. The filled 
rectangle represents a GRAPE tt pulse selective for Ci or C2 with duration of 2 ms. The phases 
of the pulses are denoted above the rectangles. 
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FIG. 3: (Color online), (a) Experimental results for error correction (EC), decoding (DE), and free 
evolution decay (FED). For each delay time, we take five data points by repeating experiments, 
shown as • for EC, O f° r DE and A for FED. The averages are shown as " x", "+" and ~k, which can 
be fitted as 0.9828-0. 0166t-0.5380t 2 +0.0014t 3 with relative fitting error 0.73%, 0.9982-0.4361t+ 
0.1679t 2 +0.2152t 3 with relative fitting error 0.57% and 1.0056 - 0.4164t + 0.3363t 2 -0.2123i 3 with 
relative fitting error 0.45%, shown as the thick dash-dotted, solid and dashed curves, respectively. 
The ratios of the first-order decay terms in the fitted curves are calculated as 26.2 ± 0.3 for DE 
and EC, and 25.0 ± 0.3 for FED and EC, respectively. The thin dash-dotted, solid and dashed 
curves show the fitting results using the ideal data points from simulation by introducing factors 
of 0.983 ± 0.006, 0.998 ± 0.007 and 1.0098 ± 0.0064 for EC, DE and FED, respectively, (b) Results 
in the previous experiment Q], shown as the data marked by "o" and "□" for EC and DE, which 
can be fitted as 0.7895 - 0.0957t - 0.0828t 2 + 0.0370t 3 and 0.8539 - 1.1021t + 0.8696t 2 + 0.0378t 3 
with relative fitting errors 0.89% and 0.98%, respectively. The ratio of the first-order decay terms 
is 11.5 ±0.2. 
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